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ABSTRACT 

We consider the shape of an accretion disc whose outer regions are misahgned with the 
spin axis of a central black hole and calculate the steady state form of the warped disc 
in the case where the viscosity and surface densities are power laws in the distance 
from the central black hole. We discuss the shape of the resulting disc in both the 
frame of the black hole and that of the outer disc. We note that some parts of the disc 
and also any companion star maybe shadowed from the central regions by the warp. 
We compute the torque on the black hole caused by the Lense-Thirring precession and 
hence compute the alignment and precession timescales. We generalise the case with 
viscosity and hence surface density independent of radius to more realistic density 
distributions for which the surface density is a decreasing function of radius. We find 
that the alignment timescale does not change greatly but the precession timescale is 
more sensitive. We also determine the effect on this timescale if we truncate the disc. 
For a given truncation radius, the the timescales are less affected for more sharply 
falling density distributions. 

Key words: accretion, accretion discs - X-rays: binaries - galaxies: active - galaxies: 
jets - quasars: general 



1 INTRODUCTION 

Observations indicate that accretion discs around black 
, holes can be warped. Warped discs have been observed in 
■ active galactic nuclei (AGN) by wate r masor observations 
' in NGC 4258 llHerrnstein. Greenhill fc Mor an 1996) and in 
the Circinus galaxy ( Greenhill et al.l (20031 ). A warped in- 
ner accretion disc might explain why radio jets from AGN 
are not perpendicular to the plane of the Galactic disc 
l|Kinnev et al.ll2000l : ISchmitt et al.ll2002l ). 

The two X-ray binaries GRO J 1655-40 and SAX J1819- 
2525 have also been observed to have jets misaligned with 
their orbital planes. For example, GRO J 1655-40 appears 
to have a binary orbit at 70° ( Greene, Bailyn & Orosz 2001;) 
and jet inclination at 85° (iHiellming fc Rupenlll995l ). This 
implies that there is a misalignment of at least 15° between 
the inclination of the black hole and the outer parts of the 
accretion disc. 

We consider a system with an accretion disc around a 
spinning black hole. The black hole spin is misaligned with 
the outer parts of the disc which we assume to be fixed 
by the plane of a binary companion. Lense-Thirring preces- 
sion drives a warp in the disc which reaches a steady state 
when the inner parts are ali gned with the black hole by the 
iBardeen fc PettersonI l| 19751) effect. 

IScheuer fc Feilerl (| 19961 ) calculated the shape of the 



steady disc and the timescale for the black hole to align on 
the assumption that the warping is gradual and that the vis- 
cosities and the surface density of the disc are independent of 
radius. Because in their analysis the orientation of the outer 
disc is fixed, the torque between the disc and the hole makes 
the blac k hole precess an d makes its spin aUgn wi th that of 
the disc (|King et al.ll2005l ). IScheuer fc r'eileil |l996l ) find that 
the alignment timescale and precession timescale are about 
the same in this case a nd this has been illustrat ed with nu- 
merical simulations bv iLodato fc Pringld l|2006l ). However, 
in more realistic discs the surface density is a decreasing 
function of radius. For idealised discs in which in which shear 
viscosity v\ varies as a power law, ui oc , in radius R the 
steady state surface density E obeys E oc R~ ^ because viT i 
tends to a constant far from the inner edge (|Pringlelll98TI ). 
Typically the power 13 lies in the range Q ^ (3 ^ 2. 

iNataraian fc Armitagel (|l999h discuss how the align- 
ment timescale depends on /3 and suggest on dimensional 
grounds that the timescale increases by about a factor of 10 
as (3 changes from to 1.5. We comp ute the shape of the disc 
but our analysis differs from that of INataraian fc Armitagg 
(1999) in that we compute how the precesssion and align- 
ment timescales vary with the density distribution {(3) when 
the warp radius and the accretion rate are fixed. We find that 
the alignment timescale does not vary strongly with (3 but 
that the precession timescale is more sensitive.. 



2 R. G. Martin, J. E. Pringle and C. A. Tout 



2 STEADY STATE SOLUTION 



Our analysis follows that of IScheuer fc Fdtol l| 19961 ). We 
consider the disc to be made up of annuli of width dR and 
mass 2nTiRdR at radius R from the central star of mass M 
with surface density E(_R, t) at time t and with angular mo- 
mentum L — {GMRY^^TjI = LI. The unit vector describing 
the direction of the angular momentum of a disc annulus is 
given by / = (IxJvJz) with \l\ = 1. 

Like IScheuer fc Feilerl (|l996l l we derive a solution for 
the steady state disc profile describ ed hy W = Ix + ily ■ We 
use equation (2.8) of lPringld (|l992l ') setting dL/dt = and 
adding a term to describe the Lense-Thirring precession to 
give 



= 
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There are two viscosities, vi corresponds to the azimuthal 
shear (the viscosity normally associated with accretion 
discs) and 1^2 corresponds to the vertical shear in the disc 
which smoothes out the twist. The second viscosity acts 
when the disc is non-planar. The Lense-Thirring precession 
is given by 



2GJ 



(2) 



IIKumar fc Pringle|[l985h where the angular momentum of 
the black hole J = Jj with j = {jx,jy,jz) and |j| = 1 can 
be expressed in terms of the dimensionless spin parameter 
a such that 



J = acM 



m 



(3) 



We take both viscosities to have power law form so that 
i^i = J'lO -5- and 1^2 = 1^20 [ ] , (4) 



where i^io, 1^20, P and 7 are all constants and -Ro is some 
fixed radius. The surface density is 



(5) 



We take the scalar product of equation ([T]) with I and find 

(6) 



RdR 



because l.dl/dR = when \l\ — 1. We assume the warp 
is gradual enough that we can neglect the non-linear term 
l.d^l/dR? = — \dl/dR\^ . We consider the effects of neglect- 
ing this term at the end of this section. This equation has 
the solution 



CoR^^^ + Ci 



(7) 



where Co and Ci are constants. We set L = at i? = and 
L = {GMR)iT, so Ci = and 

-/3 

(8) 



L = iGMR)^T.o ( -|- 



in the steady state as in the flat case. 



Substituting equation (|8]) into equation ([1} we find 
dR 



i?2 



1r r 



(9) 



We work in the frame of the black hole where J /J — (0, 0, 1) 
so that ujp = (0,0, ojp) and cjp xl = {—ujply,u!plx,0)- 

We add the a;-component of equation ([9]) to i times the 
y-component, where i = \/— 1 and set W = l^ + ily to obtain 



^R-i-^W = 4; 

dR 
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where 



so that 
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To simplify the analysis we let /3 = 7 so that the two vis- 
cosities obey the same power law and thus the ratio vi/u2 is 
independent of radius (c.f. Lodato fc Pringle, 2007). Setting 
X = we find 



3 + 2g 



(i + pr i±£ d 

—x^+iJ — 
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(13) 



and setting W{R) = R'iV{R) = ^^(i+f) v{x) we obtain 
,d^V dV 



dx"^ ^ dx 



+ 



2k 2 



(2(l + /3))2 ■ \ l + p 
This is a modified Bessel equation with solution 
2 1 \ / 2 



^^^ = 0. (14) 



V ^AI 



2(l+/3) 
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t^x- +BK, 



2(1+0) V 1 -I- /3 



(15) 

where Iu{z) and K^{z) are the modified Bessel functions 
of the first and second kind respectively and A and B axe 
constants to be determined. 

We know that ^ as i? — > because the inner disc 
is aligned with the spin of the black hole. Thus we take n to 
have a positive real part and A = so that the full solution 
is 



W 



K 



2 jR-h(^+0y], (16) 



H-/3 



In order to find B we need to consider what happens as 
7? 00. At large radius the disc tilt is taken to be fixed and 
thus we let W — > Woo which is a constant. Thus as R —> 00 
we find 



W 



BR, 



1 



K2 



' 2(l + /3) 



2{1+I3)J \l-K/3, 
where we have made use of the relation 

as 2; ^ 0. Rearranging we now find 



= Wo. 



B = 



2(1+0) 



Rn 



(17) 
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(19) 



Alignment and Precession of Black Holes 3 



and 



,2(l+,3) 



K2 

1+/9 
2 



We note that if /3 = this reduces to 



W = W^ exp 



1/2- 



(20) 



(21) 



where we have made use of the identities Ki^2{z) = 
e~ \/n/(2z) and r(l/2) = a/tt. This is the solution found 
bv lScheuer fc Feilej l|l996l ). 

T he radius wh ere the warp in the disc typically occurs, 
7?warp l|Scheuer fc Fcilcr 19961), can be found by balancing 
the terms on either side of equation ((9)1. We find 



2lu- 



V20 



P r/3 



l/(l+/3) 



(22) 



In order to compare difi'erent power laws for viscosity in a 
reeisonable way we want to keep vi, U2 and E the same at 
J?warp where the torques are greatest. We therefore set 



Rn — R^ 



(23) 



^20 V2QC'^ 

and note that 1/20 now corresponds to the value of V2 at 
the radius where the disc is warped. Equation (|20p with 



-iR^ip becomes 



W 



2VK„ 



X K 



2(1 + /J) 



1 + 



(l-») 
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(24) 



Note that we choose the negative root of —i in equation p2p 
because we want the real part of the argument of the Bessel 
function to be positive. We note that in this solution we 
have a term of the form 

A-^e-'^/^'x) = e-""/'' (ker^a:) - ikei^x)) , (25) 

where ker^ and kei,y are Kelvin functions (|Watsonlll96^ . 

The second order term, \dl/dR\'^, which we choose to 
neglect, from equation ((6)1 has a magnitude that is largest 
in the disc around -Rwarp but is negligible for small incli- 
nation angles. The largest error occurs when the outer disc 
is inclined at an angle of 7r/2 to the black hole. Then the 
relative magnitude of the neglected term is 0.049 i^2o/'^io for 
/3 = 0. For /3 = 3 it grows to 0.093 !/2o/'^io. If V20 < vio the 
analysis is good for all inclinations of the outer disc. When 
this inclination is reduced to 7r/6 the relative magnitude of 
the neglected term has fallen to 0.023 1/20/1^10 for P ~ 3 and 
0.012 z/2o/;^io for 13 = 0. 



3 THE SHAPE OF THE DISC 

In Figure [T] we plot the solution W = Ix + ily for various 
values of /3 as Ix = K(VK) against ly = Q{W) with Waa = 1- 
Note that since the problem is a linear one, we may take 
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Figure 1. The steady state disc, ly against Ix as R changes in 
the case with Woo = 1- The solid hne has /3 = 0, the dotted line 
has 13 = 3/4 and the dashed line has [3 = 3. The lines begin at 



R = 0£Ltlx 



ly 



■ and the dots on the curves are where _R = 1, 



10,100 and 1000 



VKoo = 1. In this plane, the completely flat, but inclined, 
disc would be a point at = 1. As ii ^ the disc becomes 
steadily more aligned with the black hole spin and W — > 0. 
The lower the value of /3 the less the disc is twisted. The 



and the dots on the 



lines begin at R = at Ix 
curves are where R/Rv,a.ip = 1, 10,100 and 1000. 

Kelvin functions are in effect combinations of the or- 
dinary, oscillatory Bessel functions J{x) and Y{x) and the 
modified, non-oscillatory Bessel functions I{x) and K{x). 
Because of the nature of Kelvin functions and because, as 
R 0, the argument of the Bessel function tends to infinity, 
the solution W{R) circles the origin an infinite number of 
times as _R — > while at the same time approaching the ori- 
gin exponentially. Thus, as we approach the origin, the disc 
becomes ver y twisted, b ut very fiat. This explains why, as 
remarked bv lScheuer fc F cilor (1996), numerical integration 
packages tend to fail for this problem. 

The inclination of the disc relative to the black hole spin 
direction z at radius R is 



e{R) 



(26) 



The solution for W is in the frame of the black hole where 
J/ J = (0,0,1) and we have the disc angular momentum 
vector 



I = {^{W),'=s{W), sjl - \W{R)\ 



(27) 



Thus the inclination of the disc at radius R is this frame is 

e^{R) = cos-' (Vl-iW"(^)P) • (28) 

We can also find the disc inclination in the frame aligned 
with the outer disc regions. In this frame the inclination of 
the disc tends to zero as _R — > cxd. This means as R ^ 00, 
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I (0, 0, 1). In this frame we take the angular momentum 
of the black hole, J', to have j'y — ai t — and so 

J' = J (— sinr;, 0, cos r;) (29) 

where rj is the angle of inclination of the black hole spin 
to the outer disc axis and prime denotes quantities in the 
frame of the outer disc. The disc angular momentum direc- 
tion vector becomes 

l'^ = Ix cos rj ~ sin 77 

= K(W)cosr?- - \W\^smri 

l'^ — Ix sin ri + Iz cos tj 

= 5R( VK) sin ri + ^\ - \W\^ cos 7? 

= ly (30) 

then 

6»2(i?) =cos~' {^{W) sin77+ x/1 - |M/(J?)|2 cos r?) . (31) 

If we transform back to the frame of the black hole we 
would need the transformation 

Ix ~ I'x cos rj + l'^ sin 77 
Iz — —I'x sin Tj + I'z cos 77 

ly = ly (32) 

and so we see that in the frame of the black hole as i? — > cx) 
Z (sin 77, 0, cos 77) (33) 
because I' = (0, 0, 1) and so we find 

W^oo=sin77, (34) 

or more generally 

Wo. = -ij'x+ij'y). (35) 

3.1 Application to GRO J1655-40 

As a real example we consider parameters relevant to the 
source GRO J1655-40 and take 77 — 0.2618 corresponding 
to an angle of 15°. We assume that the outer disc plane 
corresponds to that of the binary and that the inner disc 
is aligned with the spin of the hole which is parallel to the 
observed jet. In Figure[2]we plot the inclination in the frame 
of the black hole and in Figure [3] of the binary. The former 
shows that the warp steepens with increasing /3. The latter 
shows that, in the frame of the binary, the inclination of the 
disc near _R = is higher than the inclination of the black 
hole. This hump, that is due to the precession of the disc 
around the hole, could shield the binary companion from the 
black hole radiation. 
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Figure 2. The inclination of the disc against i?/i?warp in the 
frame of the black hole. The solid line has /? = 0, the dotted line 
has /3 = 3/4 and the dashed line has /3 = 3. 
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Figure 3. The inclination of the disc against -R/Rwarp in the 
frame of the binary. The solid line has (3 = 0, the dotted line has 
/3 = 3/4 and the dashed line has (3 = 3. 



4 ALIGNMENT AND PRECESSION ^ork in the frame of the black hole and the torque on the 

TIMESCALES black hole is is given by 

Having obtained the shape of the disc we are now able to r u>p x L 

calculate the mutual torque between it and the black hole. ^~dt ~ I ^ 2-KRdR (36) 

Because we assume the outer disc is fixed, this enables us to „ 

compute the timescale on which the black hole spin aligns = / ujp{—ly,lx,0)C22nR~'!~'^dR, (37) 

with it and to compute the precession rate as it does so. We JRin 
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where Rin and -Rout are the inner and outer edge of the 
disc and C2 = (G'M)^^'^Eo-Ro is a constant. Adding the x- 
component to i times the y-component we obtain 

rRaut 



/■Kout 3 

= -2niujpC2 WR'^'^dR (38) 



d{Jx + iJy) 
di 



and using equation (|16|l to substitute for W we find 
d{Jx + iJy) 



dt 



2^TiuJpC2BR^^ / R~i' 



X K 



2 jR-hi^+m] dR,{W) 



where the constant B is defined by equation (|f9|l . We sim- 
plify the expression by taking dimensions out of the integral 
to obtain 



d{Jx+iJy) 47ria;pC2 „„i / 1 '<^+« 
di =^^^-^0 



/Zout 1 + 2/3 

K ^ i (z)25TT+3y dz, (40) 



2(1 + 13) 



where 



z{R) 



-K^ R 2 



|(l+/3) 



1 + /3 

so that 2in — z{Rin) and Zout = z{Rout)- 
In Appendix 1 we show that 

r" (1 — z)oo 



(41) 



/ 

^0 



K,{z)z'' dz = 2'*-^r(i(l-c + d) 



r(-(l + c + d)) (42) 



if K(c-d) < 1 and K(c + d) > -1. So with c = 1/(2(1+/?)) 
and d = (1 + 2/3)/(2(l + /?)) this integral is valid for l3 > 
— 1/2. We now have 

3 + 413 

d{J. + iJy) ^l![!^p£^5^| C^-^\ 2" WW 

dt - n 



1 + ■■ V 2A.t 

1 + 2/3 



X r 



2(1 +/3) 



(43) 



and using equation p9|l to eliminate B we get 



d{Jx+iJy) _2TTiujpC2„, ( 1 + I3 \ '1+"^ 

dt - 1+/3 ; 

In Section[3]we found that Woo ~ —{jx + ijy)- We define 



- /3 



1+2/3 
2(l+f() 



(2(l + ;3)) 



(45) 



so that 



jx +ijy 1 + /? J V K2 

Using equation ([2} we get 



1+23 



r^^dt. (46) 



d_iM±}2A . + ^)T|.,-^r.dt (47) 

+ iJy c2 



and equation (|lip for k we get 

.?a: + ^3y C 

V 1 + 2/3 

x(i + /3)T+F^ r^dt. 



+3T 



(48) 



Then using equations ([2]) and ([Sjl we find 



l+/3\T+77 faG^M 



X r,3dt. 

We can rewrite equation (|49p as 

+ %) 



Jx + y?; 



N . 1 ,.. dt 



(49) 
(50) 



= ^5 ^0 



by setting 



V20C-' 

We can then integrate to find 



c 

~f2 71^2 



l + 2f3 



(51) 



jo: + ijy =A exp 

=yl exp 
X exp 



, , 1 t 

-(-i) 2(1+0) — 



t 

4(1 +/3); T 

TT \ t 



A{l+fi))T\' ^^^^ 

where A is the value of jx + ijj, at t = 0. In Figure 2] we 
plot the evolution of jy against jx with = 1 so j^, = 1 
and jy = Q &t t = 0. The points along the lines are at times 
t = 1, 2, 3 and 4T. 

4.1 Alignment Timescale 

Thus the timescale for alignment of the black hole is 
T 



4{l+/3) 



2 / 9 \ T+n 

2(1 + 13) I ^ \ ^'^ 



'27r(G3M)i/2Eo^° 

v^of? ) r ( 1+^^ ^ 

l2{l+/5)y 



COS (^4(i!,_^)) 

Putting _Ro = -Rwarp we find 



(53) 



47r(G3M)i/2Eo 



l/2y -"-warpV-l ~ 



(2(l + /3)) 



(2(l+/3)) (4{l+/3)) 



(54) 
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Figure 4. The angular momentum of the hole, jy against jx 
evolving in time. Initially jx = 1 and jy = 0. The points along 
the lines are at times t = 1, 2, 3 and 4T. 
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Figure 5. The alignment timescale against /3 normalised by the 
alignment timescale when /3 = 0. 



If /3 = then we get 

Malign (0) 



/ acM \ " 



(55) 



^/2^TT,o \V2oG ) 

which agrees with IScheuer &: Feileil (|l996l ) who omitted a 
factor of \/2. 

We can write the timescale to align with j3 in terms of 
the timescale with /3 = so that 



talign(/3) ^ (l+/3)-T+g 
ialign(O) 



2(l+/3) 



V2 



r ( ^±2^^ cos ( 



(56) 



In Figure [5] we plot this ra tio as a function of j3. Like 
iNataraian fc Armitagd (|l999l ). we see that as /3 increases, 
the timescale of alignment increases, but only by just under 
a factor of 2 as /? varies from to 3. 



4.2 Precession Timescale 



The precession timescale is 



T 



sin (iTTTw) 

= cot 



^ align 



^4(1 + /3)/^"^"- ^^^^ 

Thus, as found bv lScheuer fc Feileij (119961 ) . when /3 = these 
two timescales are identical. In Figure [6] we plot the ratio of 
the precession timescale to alignment timescale against /3. If 
/3 = then the alignment and precession timescales are the 
same but if /3 > then the precession timescale is longer 



02. 




Figure 6. The precession timescale divided by the alignment 
timescale against /3. 



than the alignment timescale and increases with increasing 
f3. This is also apparent from Figure [l] 



5 TRUNCATION OF THE DISC 

The results in the preceding Section were obtained under the 
assumption that the outer disc radius is infinite. We actually 
only require Rout 2> i?warp. However, this may not always 
be the case. For example in binary star systems the disc is 
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truncated at the tidal radius. We consider here the effects 
of truncating the disc at finite radius. 

The mass of the disc in steady state is given by 

/■-Rout /"-Rout 

Mdisc = 27r / T.RdR= 2tv'EoR^ / R^~'^ dR 
27rSo 



R^\2-P) 



(58) 



if /3 7^ 2. ff /9 < 2 then for finite mass -Rout must be finite 
and ff /3 > 2 then i?in / 0. 

We consider the integral in equation (|40p . If we let z = 
(1 — i)y where y is real and 



y ■ 



then the integral becomes 



(59) 



rVaut 1 + 2(3 

Q= / {l-i)K^^{{l~i)ym-i)y)^i^ dy 

J Via 



2(l + /3) 



We now define 



1 + 2|3 



(1 - ((1 - i)ym - ^)y)^^ dy 

/o°°(l - i)K_^{(l - i)y){{l - i)y)m) dy 



(60) 



(61) 



where 



J/out = 



1 + 



V20 



13 \ 2 



K^"-'' (62) 



so that in units of i?o — -Rwarp 
V2 f Ro 



2/out 



1+13 \R. 



4(l+/3) 



(63) 



Thus the quantity P(j/out) encapsulates the effect of 
truncating the disc at radius i?out(j/out). In Figure [7] we 
plot the effect of truncating the disc on the timescale for 
/3 = 0, 3/4 and 3. We plot |P] against i?out/i?warp. If |P| = 1 
then truncating the disc at that radius has no effect on the 
timescale for alignment. We see that if /3 is higher, we can 
truncate the disc closer to the central black hole without af- 
fecting the timescales. If /3 = 3 we could truncate the disc at 
Pout = 4i?warp whereas if /3 = we cannot truncate the disc 
within J?out ~ 10'' -Rwarp and leave the timescales unchanged. 



6 COUNTER ALIGNMENT 

IScheuer fc Feileil l|l996l ) considered black holes almost anti- 
parallel to the discs. In this case we reverse the sign of ujp 
and obtain 



1^20 

and hence equation (|24p becomes 
2Woo (i)^rTT7?T 



(64) 



w =- 



V2 
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X K 



2UTW1 \1 + (3 
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R 
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(65) 
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Figure 7. The modulus of P against -Rout/fiwarp. The solid line 
has /3 = 0, the dotted line has /3 = 3/4 and the dashed line has 
/3 = 3. 



and equation (|50p becomes 
d{jx + ijy) 



Jx + IJy 



. 1 dt 

■ i 2(l + /3) 

r 



(66) 



so that jx and jy increase exponentially and the disc realigns 
initially on the same timescale. 



7 CONCLUSIONS 

We have derived the steady state profile of a warped ac- 
cretion disc in the case when the viscosity and the surface 
density vary as power laws in radial distance from the cen- 
tral black hole. We fin d that, compared to the analysis of 
IScheuer fc Feileil (|l996l ) where constant surface density was 
assumed, for more realistic situations in which the surface 
density is a decreasing function of radius, the timescale for 
alignment of a black hole with its accretion disc increases 
slightly while the timescale of precession is more greatly 
increa sed. For constant surface density IScheuer fc Feilen 
(|l996l l found these two timescales to be the same. For more 
realistic density distributions we find that the black hole 
precesses at a much slower rate than the rate at which it 
aligns. If this process were responsible for changing the jet 
direction in an observed source then we would predict that 
there should be little evidence of precession. 
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APPENDIX 1 



From tables of integrals (|Gradshtevn fc RvzhikI Il980l ) we 
find 



f 

Jo 



K,{z)z'' dz = 2'*"'r Q(l ~c + d) 



xr(-(l + c + d)) (67) 



valid for 5R(c - d) < 1 and K(c + d) > -1. We need the 
integral over (0, (1 — i)oo) along the contour Co in Figure|8l 
There are no singularities in Kc{z)z'^ other than a branch 
point at the origin when d is not an integer. We have 

^(1 — i)oo f-oo 

K,(z)z-'dz= / K,(z)z-'dz + Q2 + Qz 
Jo 

(68) 

where Q2 is the integral along contour C2 over z = Se** as 
5 — > 00 and Qs is the integral along contour C3 over z — ee'* 
as e with -7r/4 < e < 0. We find 



?2 = lim 

S — *oo 



-7r/4 



(69) 



so that 



IQ2I < lim 

S — »oo 



t/4 



\K,{Se'')S''+''\de 



1/2 



lim 

S — 'OO 



-,d+ 



1/2 r 

Jo 



iv/4 



' de 



(70) 

as S — > CX3 if cos 5 so that — 7r/2 < 9 < n/2. We have 
used the asymptotic expansion (|Watsonlll966l ) 



Mz)^(^y\-il + ...), 



(71) 




Figure 8. The complex contour of integration. The contour Co 
is the integral we need. We know the integral over Ci . 



as z ^ 0. We have 



rO 

= lim / i 



iBs d+1 ie(d+l) 



de 



(72) 



and using the approximation in equation p8|l we find 



e^O 




TT r( 


c) 


4 21 


— c 








\Kc(ee"')\de 



t/4 



(73) 



because ?R(c — d) < 1. Thus 

r" (1 — i)oo 



/ 

^0 



roo 

Kc{z)z'^dz= I Kc{z)z'^dz. (74) 
Jo 



